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ABSTRACT 

We use the UKIDSS Ultra-deep survey (UDS), currently the deepest panoramic near infra-red 
survey, together with deep Subaru optical imaging to measure the clustering, number counts 
and luminosity function of galaxies at z ~ 2 selected using the BzK selection technique. 
We find that both star-forming (sBzK) and passive (pBzK) galaxies, to a magnitude limit of 
Kab < 23, are strongly clustered. The passive galaxies are the most strongly clustered pop- 
ulation, with scale lengths of ro — 15.0^2 2^^^Mpc compared with ro = 6.75^o isyh^^Mpc 
for star-forming galaxies. The direct implication is that passive galaxies inhabit the most mas- 
sive dark-matter halos, and are thus identified as the progenitors of the most massive galaxies 
at the present day. In addition, the pBzKs exhibit a sharp flattening and potential turn-over in 
their number counts, in agreement with other recent studies. This plateau cannot be explained 
by the effects of incompleteness. We conclude that only very massive galaxies are undergo- 
ing passive evolution at this early epoch, consistent with the downsizing scenario for galaxy 
evolution. Assuming a purely passive evolution for the pBzKs from their median redshift to 
the present day, their luminosity function suggests that only ^ 2.5% of present day massive 
ellipticals had a pBzK as a main progenitor. 

Key words: Infrared: galaxies - Cosmology: large-scale structure - Galaxies: High Redshift 
- Galaxies: Evolution - Galaxies: Formation. 



1 INTRODUCTION 

There is growing evidence to support the view that the most 
massive objects in the Universe were the first to as sem- 
ble and comp lete t heir star form a tion jKo dama et al. 20041 : 
Thomas et al. 20o3; lOe Lucia et al] |2004| : Isundv et all llOod : 
Stott et al. l2007i): this phenomenon has become known as down- 
siz.ine Jcowie. Songaila. & Bargeill 19991) . A number of key issues 
remain unresolved however. There are indications that the build 
up of the galaxy col our bimodality occurs around 2 ~ 2 (e.g. 
ICirasuolo etal]|2007h . but the precise evolutionary path from the 
distant Universe to the present day is still undetermined. The mech- 
anism that te rminates the major episode of star-formatio n is poorly 
understood (Benson. Kamionkowski. & Hassanil |2005|) . and it is 
now clear that massive galaxies unde rgo significant size evolution 
from z ~ 2 to the present day (e.g. ICimatti et al.ll200^ and refer- 
ences therein). 
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Galaxy clustering is an important tool for investigating these 
populations, since the amplitude of clustering on large scales (> 
IMpc) can provide a measurement of the d ark matter halo mass 
jMo & Whit j 1 19961 ; ISheth & TormenI 1 19991) . In principle, there- 
fore, one can relate galaxy populations from the distant past to the 
present day by tracing the evolution of dark matter halos within the 
context of a framework for structure formation. 

A full exploration of these issues will require large spectro- 
scopic surveys of IR-selected galaxies over a representative vol- 
ume of the distant Universe. Such surveys are currently at an early 
stage, so recent work has focused on the photometric colour selec- 
tion of passive vs. star-forming galaxies at the crucial z ^ 2 epoch. 
The t wo key methods to date are the BM/BX selection dErb et aU 
l2003h , which is an extension of the Lyman-brea k dropout tech- 
nique, and the BzK technique jPaddi et alj|2004 ') which is based 
on K-h&nd selection. Recent studies have shown that the BM/BX 
technique is reasonably efficient at selecting actively star-forming 
galaxies at z ~ 2 , but largely misses the passive galaxy p opulation 
at these epochs dOuadri et al.ll2007l : iGrazian et al.ll2007l) . In con- 
trast, the BzK method appears to be the most complete of the broad- 
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band tec hniques for the selection of both star-formins and passive 
galax ies jPaddi et al.l2004lvm Dokkum et alj2006l : rGrazian et all 
Based on initial /('-band selection, in principle it does not 
suffer the same heavy biases caused by dust or the age of the stellar 
populations. 

An alternative is to use a larger set of filters over a range of 
wavelengths and infer a galaxy's redshift from comparison of the 
calculated magnitudes with a set of templates ( Cirasuolo et al .120071 
and references therein). This work makes use of such photometric 
redshifts and in principle one can derive the stellar age of a galaxy 
using such templates. However due to the uncertainties in deter- 
mining stellar ages from the template fits, and the desire for com- 
parison with the literature, the bulk of the present analysis is based 
on simple BzK selection. 

The BzK techn iqu e was first develope d using the K20 survey 
jCimatti et ar]|2002l) bv lDaddi e,V^ ( l2004h . and preferentially se- 
lects star-forming and passive galaxies in the redshift range 1.4 < 
z < 2.5 by using the B-, z- and K^-band broadband filters. Re- 
cent work has sh own that the number counts of these populations 
differ markedly jKong et al.l l2006l : iLane et al.l |2007|) . with star- 
forming galaxies being far more abun dant at all m agnitudes. The 
abundance of sBzK galaxies allowed iKong et al] (12006,) (hence- 
forth K06) to perform a detailed clustering analysis, segregated by 
limiting K-band luminosity in the range 18.5 < K^ega < 20. 
They found that the clustering of s BzKs is strongly dependent on 
K-haiid luminosity. Hav ashi et al.l ( |2007|) studied the clustering of 
sBzKs over a smaller area (180 arcmin'^) but to much greater depth 
(Kab < 23.2) and confirmed this strong luminosity dependence. 
The sBzK population therefore appear to inhabit a range of halo 
masses, and are therefore likely to be the progenitors of a wide 
range of present day galaxies. 

There have been fewer studies of passive galaxy clustering at 
high redshift, in part because of their relatively low surface den- 
sity. K06 found early evidence that pBzKs are more strongly clus- 
tered than sBzKs, although the limited depth and area of this survey 
(320arcmin^ to Kvega = 20) gives rise to significant statistical 
uncertainly, particularly given the relatively lo w surface density of 
pBzKs at these magnitudes (0.38 arcmin" '^). iBlanc et ^ j2008l) 
measured the clustering and number counts of BzK-selected galax- 
ies to the same depth as K06, but over a much wider area (0.71 
deg^ total between two fields). Their BzK clustering amplitudes 
are smaller than those of K06 but consistent at the la level. Analy- 
sis of deeper survey data over a similar area is required to confirm 
thes e results, and we p resent such an analysis in this work. 

iLane et al.l ( |2007|) (henceforth L07) used the Early Data Re- 
lease (EDR) from the UKIDSS UDS to investigate the number 
counts and overlap in colour-selected galaxies, including those se- 
lected by the BzK technique. They found that the flattening in pBzK 
number counts observed by K06 extends to become a plateau at the 
EDR depth of Kab < 22.5 and noted that if this were to con- 
tinue to turn over it could imply an absence of high-redshift passive 
galaxies at low luminosities. 

In this work we investigate the number counts and clustering 
of BzK-selected galaxies over a substantially greater depth and area 
than any previous study. Section 2 describes the data and galaxy 
number counts. In Section 3 we present an angular clustering anal- 
ysis, which is extended by de-projection to infer the real-space cor- 
relation lengths. In Section 4 we use photometric redshifts to derive 
a luminosity function for these populations, and compare with the 
present-day galaxy luminosity function. A discussion and conclu- 
sions are presented in Section 5. Throughout this paper we assume 
a flat ACDM cosmology with = 0.3 and h = 0.73. 
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Figure 1. BzK colour-colour diagram for sources in the UDS DRl. The 
sBzK and pBzK selection regions ar e marked ac cordingly. Also visible ai'e 
the passive galaxy track identified in lLane et alj t20Q7 ), and the ste l lar lo- 
cus used to match our photometric system to that of IPaddi et alj j2004h 
(see text). We show mean error-bars for the galaxies in the BzK selection 
regions. 



2 SOURCE IDENTIFICATION AND NUMBER COUNTS 

2.1 Data set and sample definitions 

The UKIRT Infrared Deep Sky Survey (UKIDSS) has been un- 
derway since spring 2005 and comp rises 5 sub- surveys covering 
a range of areas and depths (Lawrenc e et alj2007h . This survey has 
bee n made possible by the advent of the UKIRT Wide-Field Cdxa- 
era dCasali et alj2007l) . 

We base the present study on data from the deepest compo- 
nent of UKIDSS, the Ultra Deep Survey (Almaini et al. in prepara- 
tion). This aims to reach final depths of Kab = 25.0, Hab ~ 
25.4, Jab ~ 26.0 (5a, point source, 2") over an area of 0.8 
deg^. The size of the UDS field significantly reduces the ef- 
fects of cosmic variance and on this scale is the deepest near- 
infrared survey to date. For this work we use the UDS DRl re- 
lease ( Warren et al. 2007), which reaches 5a (point source) depths 
of Kab = 23.5 and Jab = 23.7. For details of the complete- 
ness estimation, im age stacking, mosaicin g and catalogue extrac- 
tion procedures see iFoucaud et al.l ( |2007|) and Almaini et al. (in 
preparation). In addition to these data, deep B,V, R,i' , z' imag- 
ing is also available from Subaru Suprimecam to limiting depths 
of Bab = 2 8.4, Vab = 27 8, R ab = 27.7, iab = 27.7 and 
z'ab = 26.7 dFurusawa et al.ll2008l) . The UDS and Subaru survey 
areas are not entirely coincident, which reduces our usable area to 
0.63 deg^ in this analysis. 

2.2 Construction of the BzK colour-colour diagram 

The origin a l Bz K selection technique was introduced in 
iDaddi et al.l ( l2004h based on photometry in the Bessel B-band, 
Gunn 2-band and ifs -bands as defined on the FORSl, FORS2 and 



The clustering and abundance of galaxies at z ^ 2 3 



Table 1. Differential Number Counts in log(N/deg-^/0.5mag) bins for sBzK 
and pBzK in the UDS DRl. 
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ISAAC instruments at the Very Large Telescope (VLT). To correct 
our colours to these filters we used the stellar track of K06, who in 
turn used published empirical stellar spectra f rom Piclsles ("1998^ 
and s ynthetic stellar spectra from iLeieune . Cuisinier, & Buser 
( Il997h . convolved with th e filter responses of the VLT instruments 
used in lDaddietaljj2004h . These derived colours provide a conve- 
nient reference for matching stars on the well-defined BzK stellar 
locus. We use the adjusted stellar locus of K06 in the following 
way. 

The stellar locus can be split into a main branch and a "knee" 
feature, clearly visible in the lower part of figure [T] The intersec- 
tion of these branches provides a fixed point in the BzK plane 
from which we can derive a quantitative adjustment. Using a least 
squares fit to each section, we derived the following offset in z — K 
and B — z colours to t his fix ed point to convert to the photometric 
system of iDaddi et all ( |2004|) : 



(z — K)Daddi = {z - K)uDS — 0.26 
[B - z)Dadd^ = {B- Z )UDS + 0.06 



(1) 

(2) 



This correction is then applied to all sources in the BzK plane. 
We henceforth refer to BzK photome try defined in this w ay. 

The standard BzK definitions of lOaddi et alj l l2004h were then 
used to construct our BzK sample. K-band sources (> 5(t) from 
the UDS were used as the primary catalogue, with optical magni- 
tudes extracted directly from the Subaru imaging data after care- 
ful matching of the Subaru and UDS astrometric frames. Aperture 
magnitudes were then extracted using a 2-arcsec diameter. Only 
sources outside the contaminated halos of saturated optical stars 
were used. 

All K-band sources with Kab < 23.5 were used for BzK 
selection unless both the B- and z'-band magnitudes were fainter 
than the 3a limits measured on those images (0.29% of sources 
outside of contaminated regions), since these cannot be constrained 
within the BzK plane. 

Overall this procedure results in the selection of 15177 sBzKs 
(21.7% of the full sample) and 742 pBzKs (1.06 % of the sample), 
11551 and 702 respectively at Kab < 23.0. 



Figure 2. Differential number counts of passive and star-forming BzK 
galaxies with apparent K-band magnitude, compared with the full-sample 
of ii'-selected galaxies. The error-bars shown are standai'd Poisson errors 
on the counts. Two sets of values are shown for pBzKs: the standard selec- 
tion of pBzKs (filled, red triangles); and the worst case sample described 
in the text (open, green triangles). Numerical values a re reproduced in table 
[T] Literature values from K06 and lBlanc et alj j2008l) are also shown (dot- 
dashed lines and dash ed lines respective ly). The plateau identified in previ- 
ous works (K06. L07. lBlanc et"al]l20 0g^ and subsequent turn-over prior to 
our magnitude limit are apparent in the pBzK number counts. 




Figure 3. (z — K)ab colour versus K-band magnitude. The dashed line 
is our completeness limit of Kab = 23.5 for the full sample, while the 
solid line shows how the completeness limit in z affects {z — K) colour. 
We begin to be incomplete above {z — K) > 2.94. The impact of this 
incompleteness is discussed in the text. 
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2.3 Number counts 

The differential K-band number counts are illustrated in Figure [21 
and tabulated in Table 1, for both the full sample of galaxies and 
those selected as passive and star-forming BzKs. We find a steep 
rise in the counts of star-forming sBzKs toward fainter magnitudes. 
In contrast, pBzK number counts exhibit an apparent flattening at 
K ~ 21, consistent with the findings of K06 and L07. Since these 
galaxies are sampled over a relatively narrow redshift range, this 
may imply that the passive population consists largely of luminous 
galaxies at this early epoch (see Section 4). 

To investigate the reality of the turn-over, we note that the 
galaxies around and beyond this feature have very faint B-band 
magnitudes, with a substantial fraction (62%) below the Scr de- 
tection limit in this band (Bu^-a — 28.4). A non-detection in the 
_B-band does not affect their classification, however, and merely 
pushes them red-ward in B-z colour in the BzK diagram illustrated 
in Figure 1. The more worrying class are the objects which are 
below the 3a limit in z' and B, which cannot be assigned BzK 
classification (164 objects). Figure[3]shows that even within our K- 
band limit, we become incomplete above {z — K) ~ 3. We study 
the effect of this incompleteness by considering the extreme case 
in which all of the objects with B- and z'-band magnitudes fainter 
than the 3a limits (28.4 and 26.7 respectively) are pBzKs. In this 
'worst case' sample the number counts still exhibit a plateau, as 
before, but with the absence of a turn-over. 

At our conservative limit of Kab < 23 our completeness is 
> 95% for compact galaxies (Almaini et al., in prep.), so the fea- 
ture is unlikely to be due to _K'-band incompleteness unless there 
is unprecedented size-evolution in the pBzK population towards 
fainter magnitudes compared to the general population. The ex- 
tre me compactness of passive galaxies at these redshifts observed 
bv lCimatti et alj feoOSh suggests that of the two classes, the sBzKs 
should suffer more from such incompleteness. There is no evidence 
from their number counts to suggest that they are adversely affected 
in this way, so we expect that the pBzKs are likewise unaffected be- 
low this magnitude. 

Photometric errors are another source of concern. Adjusting 
the boundaries for pBzK selection by the mean photometric errors 
for sources close to the boundaries (0.1 magnitudes), we found that 
such errors had no noticeable impact on the decline in the pBzK 
population. Additionally, since the pBzK region of the colour- 
colour plane is sparsely populated, photometric errors are likely to 
scatter fainter objects into pBzK selection, rather than the opposite. 
The faint pBzK counts are thus likely to be slightly overestimated 
because of photometric errors. 

We conclude that the flattening in the pBzK population is 
likely to be real and the subsequent turn-over is probable but not 
certain. We note that the turnover corresponds to absolute magni- 
tude Mk = —23.6 at z = 1.4, which is close to the value of L* 
determined for the K- band luminosity function at these redshifts 
jCirasuolo et"d]|2007h . Taken at face value, these results are sug- 
gestive of a sharp decline in the number density of pBzK galax- 
ies towards lower luminosities, consiste nt with expec t ation from 
downsizing. As correctly pointe d out in lGrazian et alj l l2007h and 
predicted bv lDaddi et al.l ( |2004|) . this is not necessarily equivalent 
to a true decline in the number density of passive galaxies, since 
the efficiency and completeness of the pBzK technique is largely 
untested at such faint magnitudes. In particular, there is some ev- 
idence that passive galaxies can also be found among the redder 
sBzK galaxies (at large v alues of z — K), with a possible incom- 
pleteness as high as 34% iGrazian et all ( |2007|) . We find that ran- 
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Figure 4. The angular coiTelation function for our BzK-selected galaxy 
samples. The best-fitting power laws ai'e shown, with slopes fixed to the 
fiducial value of 5 = 0.8. The pBzK galaxies are very strongly clustered, 
much more so than sBzKs, indicating that they occupy the most massive 
dark matter halos at their epoch. 



domly adding an additional 34% of galaxies to the pBzK sample 
from this part of the diagram is indeed sufficient to remove the ap- 
parent turn-over, since these are predominantly the fainter sBzK 
galaxies. 

The sBzK number cou nts a re significa n tly hig her than those 
of K06, iBlanc et alj ( |2008|) and llmai et al.l ( 120081) (not plotted); 
while the pBzK number counts are lower. The UDS field is more 
than 6 times larger than that of the K 20 survey and 4 times larger 
than the one used in llmai eta i]l l2008l) . The combined field used in 
Blanc etal. (200l) is of a similar size to the UDS, cosmic variance 
is the most likely explanation for the difference in number counts 
in this case. 



3 CLUSTERING PROPERTIES 
3.1 Angular clustering 

The 2-point angular correlation function, w{9) is defined by the 
joint probability of finding two g alaxies in solid angle elements 
50,1 and 50,2 at a given separation ( |Peebles|[l98Cll) . 



5P = n''50i502{l + w{ei2)). 



(3) 



To estimate the co rrelation function we use the estimator of 
lLandv&Szaiavl ( ll993l) . 



DD - 2DR + RR 
RR 



(4) 



where DD, DR and RR are the counts of data-data, data- 
random and random-random pairs respectively at angular separa- 
tion 6, normalised by the total number possible. Although this esti- 
mator is relatively robust against systematic errors, there remains 
a small bias due to the finite field size, which is corrected for 
an integral coristraint by a constant, C. We follow the method of 
lRoche&Eale3 ( ll999l) by using the random-random counts to esti- 
mate the size of this bias: 



C 



ENRR{e)e-° 



(5) 
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Figure 5. Photometric redshift distributions for passive (left) and star- 
forming (riglit) BzK-selected galaxies (solid line histogram). The over- 
plotted dashed lines are distributions with photometric errors de-convolved 
(see text). 



where the sums extend to the largest separations within the 

field. 

The sBzK and pBzK samples were selected to a limit of 
Kab < 23. We adopted this conservative magnitude limit to en- 
sure the minimum contamination by spurious sources. We fit a sin- 
gle power law of the form w{9) = A{9~^ — C) to the corrected 
data over the separation range 0.01 — 0.1 degrees, fixing the slope 
to the fiducial 5 = 0.8 and minimising the x^- The errors on the 
measurements are a combination of those due to shot noise (esti- 
mated by bootstrap resampling) and an estimate of those due to 
cosmic variance. The estimates for cosmic variance were found by 
splitting the field into 4 and computing the variance of w{9). 

In order to fit the clustering reliably, we wish to avoid small- 
scale excesses due to multiple galaxy occupation of a single halo. 
The lower bound was chosen to correspond to ~ 0.9Mpc (co- 
moving) at z~ 2 for this reason. The upper bound is a conservative 
estimate of the limit to which our data have enough signal for a 
reliable fit to be obtained. 

For sBzK galaxies the amplitude was found to be A — 
(deg.'''*)and for the pBzK population A — 



l-7Ql:o.ll X 10 
6.37ll;54 X 10 



^ (deg."'*). The clustering amplitude of the pBzKs 
is inconsistent with the clustering of sBzK galaxies at the 3-sigma 
level. Figure|4]shows the clustering measurements corrected for the 
integral constraint for the sBzK and pBzK galaxies. 



3.2 De-projected clustering amplitude 

The real space clustering and projected clus tering are linked by the 
relativistic Limber equation ( (Limbei|[T954 ^. If the redshift distri- 
bution of a sample is known, the Limber equation can be inverted 
and the correlation length, rp, can be calculated in a robust manner 
jPeeblesll 1 9S(i iMagliocchetti & Maddo?3l 1 999h 

To estimate the redshift distribution we use photometric red- 
shifts based on the Subaru and UDS bands previously introduced, 
with th e addition of Spitzer data taken as part of the SWIRE 
survey ([Lonsdale et al. I l2003h . The method is described fully in 
ICirasuolo et al. (2007), and consists of minimising the of syn- 
thetic galaxy templates. The distributions are shown in figure [5] 
The photometric redshift distribution was then used directly in 
the inverted Limber's equation during the calculation of tq. The 
values obtained in this manner are ro = Ib.Stlpr^ Mpc and 
7.69toAih'^ Mpc for A" < 23 pBzKs and sBzKs respectively. 
The quoted errors are due to the error in the fit to the clustering am- 
plitude and therefore take into account the shot noise and cosmic 
variance. 

The photometric redshifts are subject to errors (a/{l + z) = 
0.095 and 0.105 for pBzKs and sBzKs respectively), however, and 
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Figure 6. The dependence of clustering strength of sBzK-selected galax- 
ies on limiting Kab -band magnitude. Our meas urements (open sq uares) 
are shown together with literature values from iBlanc et all j20O8i) and 
iHavashi et aD j2007h (open triangle and asterisks respectively). Our values 
confirm the magnitude dependence of sBzK clustering strength. 



Table 2. Values of ro for sBzKs in the UDS by limiting K^g-band mag- 
nitude along with the number of objects brighter than the relevant limit. As 
for the full samples, the error ranges are derived from estimates of the shot 
noise and cosmic variance. 
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assuming that the true redshift distribution is highly peaked, these 
errors are likely to broaden the measured distribution. A broader 
distribution will result in a larger inferred clustering scale length. It 
is therefore important that we take such errors into account. We do 
so by assuming the errors are Gaussian, and then deconvolve the 
errors from the redshift distribution using a Fourier-based Wiener 
filter. A more detailed description of this method is provided in 
the appendix. The resulting de-convolved redshift distributions are 
shown in figure |5] 

The scaling lengths recovered using the corrected redshift dis- 
tribution are as follows: 15.012:2'*"^ Mpc and G.rbto.lrh'^ Mpc 
for pBzKs and sBzKs respectively. 

However, we note that our photometric redshift code is largely 
untested for pBzK and sBzK galaxies. This work will be refined by 
the use of an ongoing ESO Large Programme using VIMOS and 
F0RS2 on the VLT to target one sixth of the UDS DRl galaxies 
with photometric redshifts > 1. 

We also confirm previous cla ims for a dependence of sBzK 
clustering on apparent magnitude dHavash i et al. 2007). Figure |6] 
and Table|2]show the values for sBzKs with varying limiting mag- 
nitude. The dependence is much stronger at magnitudes of K < 22, 
indicating a strong correlation between halo mass and sBzK lumi- 
nosity for these objects. 
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3.3 BzK galaxies selected atl.4 < 2 < 2.5 

The BzK selection was defined by IPaddi et alj ( |2004|) to isolate 
galaxies in the redshi ft range 1.4 < z < 2.5. Clearly from figure 
[S] and as expected bv lPaddi et alj ( |2004|) . there are contaminating 
objects from outside of this range. Using our photometric redshifts 
we can asses how successful the BzK selection technique is in re- 
producing the clustering of objects within the desired range. 

Following the same method outlined for the full samples, 
we compute clustering amplitudes find correlation lengths of 
ll.lt\lh-^ Mpc and ^A&to^^h^^ Mpc for lA < z < 2.5 
pBzKs and sBzKs respectively. These values are slightly lower than 
those for the full samples and indicate that the high redshift tails in 
the full sample are at least as highly clustered as those within the 

1.4 < z < 2.5 range. However, the conclusions drawn from the 
full sample are still valid, namely that the pBzK galaxies are sig- 
nificantly more strongly clustered than the sBzK galaxies. 



4 LUMINOSITY FUNCTION 

A luminosity function was construc ted for our BzK-selecte d galax- 
ies in the same way as detailed in ICirasuolo et al.l ( l2007h . by us- 
ing the l/Vmax method ( ISchmidllll968h . FigurelTlshows the lumi- 
nosity function for the BzK galaxies with photometric redshifts in 
the range 1.4 < z < 2.5 (points with error-bars), compared with 
all K-selected galaxies in the same range (solid line). Al s o plot - 
ted is the 2 = luminosity function from lKochanek et"ail j200lh . 
It is clear that sBzK galaxies sample a wide range in luminosity, 
while the pBzK population are dominated by bright galaxies with 
Mk > —23. At this epoch, however, the bright end of the lumi- 
nosity fu nction is still dominate d by star-forming objects, consis- 
tent with lCirasuolo et ^ | |2007|) who found that the galaxy colour 
bimodality is weak at these redshifts. The pBzK galaxies are nev- 
ertheless likely to be among the most massive systems, since the 
mass-to-light ratio will be significantly lower for the actively star- 
forming systems. 

Under the assumption that the pBzKs passively evolve to the 
present day (with minimal merging) one can estimate their contri- 
bution to the bright end of the present-day luminosity function. The 
implied evolution is modelled from z ~ 1.60 (the median value 
for pBzKs) to z = as suming the spectral evolution models of 
iBruzual & Charloll j2003h . For simplicity we chose a solar metalic- 
ity model and a Salpeter initial mass function (Salpeter 1955) and 
assume only passive evolution from initial formation bursts at red- 
shifts 2/ = 3 and Zf = 10. This results in a Johnson A'-band 
absolute magnitude evolution in the range 0.96 < AMk < 1.21. 
Taking the mean value we can estimate the minimum number of 
present day bright {K < —23.4) galaxies that could previously 
have been pBzKs. Under our assumptions we find that 2.5% of 
such galaxies can be explained by passively evolved pBzKs. It 
is assumed that the remainder is made up of the descendants of 
bright sBzKs and merger remnants from within and between the 
two classes. 



5 DISCUSSION AND CONCLUSIONS 

We present the number counts, clustering and luminosity function 
of galaxies at 2 ~ 2 selected using the BzK selection criteria. The 
pBzK galaxies show a marked flattening in their number counts 
which cannot be explained by the effects of incompleteness, and 
a possible turn over at faint magnitudes. We conclude that, at this 



epoch, it is generally luminous, massive galaxies that are undergo- 
ing passive evolution. This is consistent with the down-sizing sce- 
nario, in which the most massive galaxies are formed first and are 
the first to evolve onto the red sequence ( Kodama et al. 2004). 

The angular clustering of the passive galaxy sample is very 
strong, approximately 4 times the amplitude of the sBzK popula- 
tion. This is in part due to their relatively narrow redshift distri- 
bution (K06), but a large difference remains after de-projection to 
the real-space correlation length. We find ro = IS.Olj'j'^^^Mpc 
and ro = 6.75to;37ft"^Mpc for the pBzK and sBzK galaxies 
respectively. Our value for the correlation le ngth of sBzKs is al- 
most twice that found bv lHavashi et al.l l l2007h . who used a sample 
with similar limiting magnitude but with smaller field size (180 
arcmin'^, compared with ~ 2250 arcmin^ for the UDS). Our sam- 
pl e consists of more tha n 10 times the number of sBzKs than that 
of iHavashi et"d] ilOOf ) and in addition we make fewer assump- 
tions regarding their redshift distribution. Further surveys reaching 
depths of KAB.iim ^ 23 are required to fully account the effects 
of cosmic variance, however. 

We also confirm that the scale-length for the clustering of 
sBzK galaxies is dependent on apparent magnitude, consistent with 
the work of Havashi et al. (2007) (figure [6]l. This dependence is 
far more significant at magnitudes below Kab ~ 22 indicating a 
strong correlation between such galaxies and the mass of the host- 
ing halo. 

In addition to this luminosity dependence there is a clear en- 
hancement at small scales in the sBzK clustering. This enhance- 
ment is indicative of multiple sBzK galaxies occupying a single 
dark matter halo. The scale at which this turn-off occurs can pro- 
vide an indication of the size, and hence mass, of the hosting dark 
matter halos. The turn-off occurs at '--^ 0.01 deg, which corresponds 
to 0.3 Mpc at z= 1.65. Halos of this size have masses between 
10" and 1O^^M0 ( IMo & Whitell2002l) . A full consideration of 
this enhancement within the halo occupation distribution frame- 
work will be presented in a future paper. The host halo mass can 
also be derived from comparing the de-projected clustering scale 
length with mod els of dark ma t ter hal o clustering evolution. Using 
models based on IMo & Whitel j2002l) we find a typical halo mass 
of ~ 6 X 10" Mq. 

Applying the model of clustering evolution to the pBzKs we 
can qualitatively conclude that they are inhabitants of the most mas- 
sive halos at their epoch (in excess of IO^^Mq); halos which will 
eventually become massive groups and clusters by the present day. 
The evolutionary path that pBzKs take with redshift will also be the 
subject of future work, using further colour selection techniques. 

Our conclusion based on the ro measurement is strengthened 
by their luminosity function. Even under the strict assumption of 
purely passive evolution, the descendants of pBzKs occupy the 
bright end of the luminosity function. Such galaxies are group 
and cluster dominant galaxies in the local universe. The brightest 
{K < 21) sBzKs have clustering scale lengths comparable to, or 
greater than, that which we have found for the pBzKs. This find- 
ing indicates that such galaxies will also become group and cluster 
dominant galaxies by 2 = 0. The implication is that we are indeed 
witnessing the epoch at which the build up of the red sequence be- 
gins. 
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Figure 7. Left: Luminosity functions for our pBzK and sBzK samples (filled and open circles respectively), plus all it-selected 
galaxies to a limit of A" < 23 (restricted to galaxies with photometric redshifts in the range 1.4 < z < 2.5, solid line). The sBzKs 
are representative of the overall population in the range 1.4 < 2 < 2.5, while the pBzKs are found to be exclusively bright objects. 
Also shown is the z = luminosity function for K-selected galaxies from[Kochanek et al. (2001) (dashed Hne). Right: The solid line 
shows the same luminosity function for 2 = gala xies as the da shed line in the right hand panel. The dashed curve shows the 2 = 
luminosity function for early type galaxies from Kochanek et aP i2001i) . while the points are the 2 = luminosity function for pBzKs 
under the assumption of passive evolution (see text). 
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filter tends to (fc) as in the simple deconvolution above. SNR is 
estimated by fitting the power spectrum of the redshift distribution 
to a function of the form a x lO"* * + c, with c identified as being 
the noise level. The resultant redshift distributions are shown in 
figure|5] 



APPENDIX 

In this Appendix we describe the Fourier method used to decon- 
volve photometric redshift errors from the n(z) distribution, which 
can then be used to invert the Limber equation. The errors in the 
photometric redshift distribution are assumed to be Gaussian and it 
is also assumed that the measured distribution is simply the convo- 
lution of these errors with the true redshift distribution. Under these 
assumptions it should be possible to deconvolve the errors from the 
measured distribution, using the fact that a convolution is simply 
a multiplication in the Fourier domain. Dividing the Fourier trans- 
form of the measured distribution by that of the Gaussian errors, 
should then give us an estimate of the true distribution. 

In practice, a relatively large number of terms in the discreet 
Fourier transform are required to reproduce the redshift distribution 
accurately, and as the the Fourier transform of a Gaussian is also a 
Gaussian, small levels of high frequency noise can be amplified 
greatly to give a spurious result. One way to avoid such a problem 
is to use a Wiener filter: 



W{k) = 



H{k) \H{kY + ^ 



(6) 



Where H(k) is the Fourier transform of the Gaussian errors and 
SNR the signal to noise ratio. In the limit of SNR being infinite, this 
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